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Abstract 

The low energy effective field model for the multilayer graphene (at ABC stacking) is considered. We 
calculate the effective action in the presence of constant external magnetic field B (normal to the graphene 
sheet) . We also calculate the first two corrections to this effective action caused by the in-plane electric field 
E at E/B <C 1 and discuss the magnetoelectric effect. In addition, we calculate the imaginary part of the 
effective action in the presence of constant electric field E and the lowest order correction to it due to the 
magnetic field (B/E < I). 



1. Introduction 

Physics of graphene demonstrates numerous deep relations with fundamental physics such as relativistic 
quantum mechanics and quantum field theory [![. The charge carriers in single- layer graphene are topo- 
logically protected gapless fermions, which in the vicinity of the nodes automatically acquire the properties 
of massless Dirac fermions with relativistic spectrum. This naturally induces Lorentz boosts, e.g., in the 
problem of graphene in crossed electric and magnetic fields @, ||. Bilayer and rhombohedral (ABC- 
stacking) trilayer [5] graphene demonstrate more exotic non-Lorentz-invariant physics of the topologically 
protected massless chiral fermions with quadractic and cubic dispersion laws, respectively. Such theories 
are intensively studied now, in particular, in a context of Hofava quantum gravity with anisotropic scaling 
(this analogy was recently discussed in Ref. 0j)- Quantum electrodynamics of bilayer graphene, which 
experiences the phenomenon of anisotropic scaling, has been considered in Ref. Q based on a semiclassical 
approach. Here we derive the effective action of electromagnetic field in bilayer and multilayer graphene 
based on exact solution of Schrodinger equation. For the case of purely electric field, this has been done 
recently in Ref. Q. We will generalize this consideration to the case of crossed fields. This allows us to 
discuss the magnetoelectric effect in graphene, that is, a change of magnetization due to electric field and 
change of dielectric polarization due to magnetic field. 

In Hofava-Lifshitz quantum gravity Q, the vacuum state is characterized by anisotropic scaling, i.e. the 
action for gravity is invariant under transformation r — > for, t — ¥ b z t, where integer z ^ 1 is the analog of the 
dynamical critical exponent in the theory of phase transitions. The Hofava-Lifshitz gravity in 3+1 spacetime 
is asymptotic safe, i.e. does not suffer from the ultraviolet divergences, if z = 3. It is instructive to extend 
the investigation of the consequences of the anisotropic scaling to quantum field theories with general z. It 
appears that such theories are relevant for the models of the multilayer graphene. These models have nodes 
in the fermionic spectrum, which are protected by the integer momentum-space topological invariant N 
expressed in terms of the Hamiltonian, see e.g. [t| (or more generally in terms of the Green's function [lojV 
Close to such a node fermions behave as 2+1 massless Dirac particles with energy spectrum £{p) = ±vp N , 
which obeys the anisotropic scaling with z — N. That is why one may expect that the effective action for 
quantum field theories emerging in such systems contains the terms obeying this scaling law, and we consider 
here such terms in the effective Eulcr-Heisenberg action for the quantum electrodynamics of the multilayer 
graphene. 
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In the simple model with J ABC-stacked graphene layers considered in Ref. Q , the topological invariant 
coincides with the number of layers, N = J. We keep in mind such model, though our consideration can be 
easily extended to the more general cases with N ^ J. We considered the Euler-Heisenberg effective action 
5 e //(B,E) in two limiting cases, electric field dominated B << E and magnetic field dominated E « B. 
Our results for the effective electromagnetic action with E « B in the system with number of graphene 
layers in the range 1 < J < 10 are accumulated in Table Q] This Table demonstrates some peculiar properties 
of the action. In particular, the ~ E 2 correction due to electric field decreases quickly with the increase of 
the number of layers while the coefficient before the main magnetic term increases. At some special value 
of J (at J = 4) the ~ E 2 correction to the effective action vanishes. There are also some special values of 
J, at which the logarithmically divergent term appears; these are J = 2, 6, 10. This situation is very similar 
to what occurs in the relativistic theories, where the logarithmically divergent terms appear only for some 



special values of space dimension D (see e.g. This suggests specific regularities in the behavior of 



quantum field theories, if they are extended to general values of D and z. 

Table Q] demonstrates that the linear term in electric field arises in the action. It is related with the 
degeneracy of the first J Landau levels at E = (they are all zero modes) [1? ]. It describes the linear 
Stark effect (similar to that in hydrogen atom, due to degeneracy of the states with different orbital quantum 
numbers), which reflects the spontaneous (broken symmetry) electric polarization emerging in the system. 
This term can be presented as a scalar product E • P, where the vector P is directed along the spontaneous 
polarization, which in the presence of electric field is oriented along the field. This term is responsible for 
the magnetoelectric effect in Eq. (|59p. 

In the electric field dominated regime B « E, our main result is the expression for the imaginary part 
of the effective action and the corresponding Schwinger pair creation rate given by Eqs. {55j) and ([70]) . 

The paper is organized as follows. In Section 2 we consider the system in pure magnetic field. In Section 
3 the system is considered in the presence of both electric and magnetic fields with E « B. In Section 4 
the opposite case B << E is considered. In Section 5 we end with the conclusions. 



2. The system in pure magnetic field 

The Euler-Heisenberg effective action S e ff (and Lagrangian L e tf) in the considered effective field theory 
of multilayer graphene is given by the expression for the partition function of the system in the presence of 
magnetic field B and electric field E: 

Z{T) = J ##exp(i J dtd 2 xi>\id t -H[E,B,t}}^ = e iTL -ss\E,B\ = e iS eff [E,B] ^ 

Here i/j is the fermion field (with hidden spin index A = 1, 2 and the index a = 1, 2 that enumerates the 
number of 2 - spinors), T-L is the one - particle Hamiltonian in the presence of external electric and magnetic 
fields E and B, T — > oo is time extent. The boundary conditions relate the values of fields at t — and t = T. 
Usually, for the solution of stationary problemsjthe anti - periodic boundary conditions are implied for the 
fermion systems (see, for example, Refs.[l2| or [HI]). However, in the case of constant nonzero electric field 
the system is not stationary. We overcome this difficulty considering the path integral in moving reference 
frame, where electric field vanishes. In this reference frame the problem becomes stationary, and we apply 
the anti - periodic boundary conditions for the fermion fields (for the details see Appendix A). 

2.1. Schrddinger equation 

We start with the case of purely magnetic field. First, let us consider the one-particle problem. Its 
solution is well known (see Chapter 2 in the book [l[ and references therein) but it is convenient to repeat 
briefly the main results for the case of J-layer graphene-like material with arbitrary J. Here we consider an 
ideal system, in which the fermionic spectrum is characterized by the nonlinear touching point £(p) = ±vp J 
(see Fig. [T] for J = 3). It roughly corresponds to the case of rhombohedral stacking of graphene layers 
in approximation, in which only the largest in-plane and out-of-plane hopping parameters are taken into 
account and the trigonal warping is neglected [lj. The spectrum with nonlinear touching points can be 
realized also in artificial materials as it is done already for J = 1 3, EH. It can also arise in relativistic 



quantum field theories where the phenomenon of reentrant violation of Lorentz symmetry takes place: the 
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fundamental Lorentz violation at high energies triggers a reentrant violation of Lorentz invariance below 
some low-energy scale leading to the non-linear touching point (see, e.g., Sec. 12.4 in the book [HI). The 
Dirac vacuum of fermions, whose massless branch of spectrum has non-linear touching points £{p) = ±vp J , 
induce the non-analytic terms in quantum electrodynamics, which obey the anisotropic scaling r — > ar, 
t -+ a J t, i.e. the magnetic and electric fields scale as B — ¥ Ba~ 2 and E -+ Ea~( J+1 h Here we will be mainly 
interested in these non-analytic terms, while the massive Dirac fields give rise to the conventional analytic 
terms in the Euler-Heisenberg [lj} action of the type ~ _B 4 /m 5 , where m is the mass of Dirac particles on 
branches with energy gap (see Eq. (46) of ll|). 



N = 1 + 1 + 1 =3 

AM-1 N=+l N=+l 




E=v F p E=v F p E=vpp 



Figure 1: Illustration of branches of spectrum emerging in the system simulating the trilayer graphene with rhombohedral 
stacking according to Ref. Q. left: Noninteracting layers. Each layer has its own Dirac fermion with linear spectrum 
£(p) = ivpp, (on this Figure vp is denoted by c, only one conical point at p = is considered in each layer). The conical 
point is characterized by the momentum-space topological charge 2V = 1, so that the total charge of the points at p = is 
N = 3 (for Dirac points in the other valley the charge is N = —3). right: Interacting layers. The spectrum is shown as function 
of p x at p y = 0. Four branches of spectrum correspond to two massive Dirac fermions whose topological charge N = 0, while 
two branches with the lowest energy correspond to the massless Dirac particle with nonlinear dispersion £(p) = ivp 3 , which 
is characterized by topological charge N = 3. For the similar stacking of J layers the branch of massless Dirac particles has 
dispersion £{p) = ivp J , which corresponds to topological charge N = J. 



We deal with the two-component spinors describing electron states in graphene placed in the external 
magnetic field B directed along z-axis (normal to the graphene plane). We consider the external vector po- 
tential in the form: A v = Bx (Landau gauge) . The one-particle Hamiltonian in a subsequent parametrization 



has the form 18 



H = v 



( 







(px + i{Py + Bx i) 



(fix -i(p y + Bx)\ 



(2) 







Here p x — —id x , v is a constant that is equal to Fermi velocity vf for the case of monolayer, l/2m (m is the 
effective mass) for the case of bilayer, etc.; we will use the units H = e = 1. This Hamiltonian describes the 
lowest energy bands at \£\ <C |tj_|, where t± is the interlayer hopping parameter and the energy £ is counted 
from the band crossing point (neutrality point). Therefore t± ps 0.4 eV plays the role of the ultraviolet (UV) 
cutoff energy (further we denote it by A) for J > 1; for J = 1 the Dirac model is applicable till the energies 
of the order of the in-plane nearest-neighbor hopping parameter t 3 eV, which gives us the cutoff energy 



for this case. The values of v and A = tj_ are related as v w A 



i-J„,J 



From the opposite side of the small 



energies, the applicability of the model @ is restricted by trigonal warping effects due to farther interlayer 



hopping terms; the corresponding energy scale is about 1 to 10 meV [l|, [18, 2l|. For the case of bilayer, it 
was shown both experimentally [21( and theoretically (see the recent paper [22j and references therein) that 
many-body effects can result into a reconstruction of the ground state at these small energies. Since our 
model is, anyway, inapplicable there we will assume an infrared cutoff (when necessary) of the order of 1 to 
10 meV. Many-body effects can also lead to a reconstruction of the ground state for the case of rhombohedral 
(ABC) trilayer, see, e.g., Ref. [23[. Detailed nature of these instabilities is still unknown and is a subject of 
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intensive debates. We will not take into account the many-body effects in this paper restricting ourselves by 
the calculations of the effective electromagnetic action in the approximation of noninteracting fcrmions only. 
Stationary Schrodinger equation has the usual form 

= (3) 

Later on we imply periodic boundary conditions in space coordinates. That is why "J" can be decomposed 
into the sum over the quantized y - momentum: \I>(i, x) = eWvV ~ l£t '4>py,£( x )- ^p v ,£ ( x ) ^ s the eigenfunction 

Py 

of the Hamiltonian (£ is the eigenvalue): 



v 



( 



[p x -i(p v + Bx) 

(px + i(Py + 



i>Py,s( x )= £ '&Pv,s( x ) ( 4 ) 



2.2. Landau levels 

Let us now introduce the notations: 

£ 



-, u = VB(x+ Py /B) (5) 



vB J / 2 : 
Then 

{-i)- J eipi = [d u + u] J ip 2 

{-%)- J &h = [du-uY^ (6) 

where = (V , i,V , 2) T - For ^1,2 we have: 

e 2 ^2^[du-u\ , [u + d u ] J 4> 2 (7) 
As for the oscillator we introduce the annihilation and creation operators a = -j%[u + d u ], a) = -^[u — d u ] 
and N = a) a. Then 

2-J/2 

e 2 i> 2 = 2 J [a)] J a J i) 2 = 2 J N(N - 1)...(JV — J + 1)^2 (8) 



The normalized solutions are 



with integer n> J and the eigenvalues e = ±y/2 J n(n — l)...(n — J + 1). The solutions for J > n > have 
tpi = and correspond to e = 0: 

V>i = 

V-2 = -^=e-» 2 /2 ffn(u) (10) 

That is why at each value of p y we have Landau levels with £ — (the degree of degeneracy is J per flux 
quantum) and the nondegenerate levels 



£ = ±vB J/2 ^2 J n(n - l)...(n — J + 1) 

As usual, we need —L x /2 < p y /B < L x /2,p y — j^Ky, K y e Z, where L XjV is the linear size of the graphene 
sheet in x, y direction. 
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2.3. Naive expression for the effective action 

For pure magnetic field at zero temperature the effective action can be written as follows (for the definition 
of the effective action S e f / and effective Lagrangian L e f f and their relation to the spectrum of one - particle 
excitations see Appendix A) 



-S efs = -L eff T =-Tj2 1^ 1/2 = -g s gv TLx 2 ^ vB v(2B) J / 2 ]T y/n(n - l)...(n - J + 1) 

n—J 



Here |£| is the sum of the energies over all energy levels; g s = 2 and g v = 2 are spin and valley degeneracies, 
T is the period of the motion. This sum is formally divergent. Using the Hawking trick of the zeta-function 
regularisation [2J| we represent it formally for J = 1 through the ^-function: 



Seff = -g s9v ^^v F (2B)^a-l/2) = gsgv l£f VFi 2B)^am 



2tt 



where we use the Riemann identity 25j, [26 1 



T (s) ( (s) cos (ttb/2) = tt s C (f - s) 



(11) 



(12) 



For J = 1 the effective action has to be Lorentz invariant. This means that constant magnetic and 
electric (E) fields must enter via the combination B 2 — E 2 /vp. Thus, we arrive at the effective action in the 
presence of both fields: 

- S eff = -L gsgv Ihh> VF [B 2 - E 2 /v 2 F f*C(Z/2) (13) 

This gives us a contribution of vacuum fluctuations (Euler-Heisenberg action 17] ) for the single-layer 
graphene. 

From here we get at B = the imaginary effective action which describes a spontaneous creation of 
electron-hole pairs from vacuum (Schwinger effect [13, M, M, M, HE M, SI H ) : 



TT T 

2ln,S eff =g s g v ±^p-E^ 2 am 

4Vp TT 2 



(14) 



This result coincides with that obtained for the vacuum persistence probability using the other methods (see 
Rcfs.|30. 7] and references therein). 

2.4- Zeta function regularization of the partition function 



Let us evaluate the Euler-Heisenberg effective action using the zeta-function regularization [351 ] . As it is 
shown in Appendix D this regularization gives the correct values of the Euler-Heisenberg effective action if 
the divergences are at most logarithmic. In the presence of the divergences that are stronger than logarithmic 
the zeta - regularized expression for the effective action gives subdominant terms. 

To use this regularization we need to calculate the determinant of a positive defined operator. Let us use 
the valley degeneracy for this purpose. We have 

(Det [id t - W^J 2 = Det (id t - Hj (a 3 [id t - H]^) = Det [id t - H][id t +U} = Det [{id t ) 2 - U 2 \ (15) 
After the Wick rotation we get 

(Det [d T - H]j 2 = Dot (-[{id T ) 2 + U 2 ]) (16) 
Zeta-function regularization (3~5| then reads: 



L. 



eff 



-rTrlog([a r -H]) 5s9 "=^9 s (^^^- 1 Tr e -[^) 2 +« 2 ] t )| s=0 (17) 
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Here the system is considered in imaginary time r € [0,7], and T = 1/T^0is the inverse temperature. 
Below we omit mentioning the exact meaning of T (real or imaginary time) . It will always be clear from the 
context. 

It is worth mentioning that the zero modes are to be omitted here. We have: 



LtJjh „ „ / 1 

eff = 



g.g v Bd.(±- I dtt^^e- [mi+1/2)f+£ - ] %=o 

v ' m,n 



AttT 



At zero temperature in the sum over m only the term with m = survives. Zero modes of % are omitted 
and we get: 



The last equation defines analytical function of s at s > s for some s . This function has to be continued 
to s = 0. 

Let us check ourselves considering the case J = 1: 

n> 1 



L X L 



lg s g v Bd s ( K ^-(2v 2 F B) 1 / 2 - s r( S - l/2)C(s - l/2))| s =o 



= h^g s g vVF B^am (19) 



which coincides with the result obtained above. 
For arbitrary J we get: 



The resulting expression for the effective action is 

- L eff = ^L 9s g v 2 J /\BV+W d s (J^ V -^B- J T(8 1/2) fj(a 1/2)) | s=0 , 
fj(s) = f^r / duu*- 1 V e -<^", a > So 



(21) 



Here the function fj(s) is defined by the given integral for s > s = 1/J while its value at s = —1/2 enters 
the expression for the effective action. The function fj(s) can be considered as the generalization of Riemann 
zeta- function. For s > 1/J it is given by the series: 



n>J 

(n)-j = n(n- l)...(n - J+ 1) (22) 
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Effective lagrangian —leff/{9s9v) 
0.04679084006 x v F [B 2 - E 2 /v 2 F f /i 



0.01989436789 x vB 2 log(^-) - 0.1125395395 x EB 1/2 - 0.03490053651 x §-, 



i D 



0.2431821436 x v 2 F B 2 A~ 1 - 0.07113483618 x vB 5/2 - 0.1949242002 x EB 1/2 - 0.02123563323 x -g™ 



0.2750863237 x v 2 F B 2 A- 1 + 0.1746037024 x vB A - 0.3462141635 x EB 



0.3484677731 x v 2 F B 2 A~ 1 + 1.486414997 x vB 7/2 - 0.4740836354 x EB 1/2 - 0.003908473257 x 



0.4434404060 x v 2 F B 2 A~ 1 - 0.3133362943 x vB 4 log( ) - 0.6561212594 x EB 1/2 - 0.001346430685 x 



vB 3 



0.5560282734 x v 2 F B 2 A~ 1 + 7.837230303^^ - 3.222245513uB 9/2 - 0.8183430244£.B 1/2 - 0.0004176771170^1^ 



0.6849063582 x v 2 F B 2 A- 1 + 7.438809597^^ + 27.35355260 x vB 5 - 1.026656793 x EB 1/2 - 0.0001187738012 x 



0.8295068879 x v 2 F B 2 A- 1 + 8.754557174^4 + 264.3063672«B 11/2 - 1.217484080£B 1/2 - 0.00003135755874- 



10 



0.9895460087^S 2 A- 1 + 10.93814698^1. - 40.348264861og( 21 ' .f )vB 6 - 1.449119539 EB 1/2 - 0.775944555 10" 5 ^ 



Table 1: Effective lagrangian —L e f j / (L x L y g s g v ). Here A is the ultraviolet cutoff (v £S A 1 J w^,), while fi is the cutoff parameter 
of the regularization scheme. Both A and ft depend on microscopic physics and should be considered as fitting parameters. 



We have (s 0): 



Us -1/2) « /j 1 ( V2) +/f(-l/2) (23) 
s 



Procedure for the calculation of the values of fj 1,0 - ) (— 1/2) is given in Appendix C. It has been shown 

that / ( 7 _1) (-l/2) may be nonzero for J = 2(2K + 1),K E Z. 

As a result for the values of J 7^ 2(2_K" + 1) we obtain L e ff = L x L y l e ff with 

-l eff = vBW^aj, 

a J = (^^.9 s .9«2 J / 2 r(-l/2)/. 7 (-l/2)), J = 1, 3, 4, 5, 7, 8, 9, 11, 12, 13, 15, ... (24) 

Now using the data of Table [2] from Appendix C one can easily calculate values of aj. 
For J = 2,6, 10, 14, ... we have (j^j « s + 7s 2 ): 

2dJ, 



«r - ^ (J/2+1) 4- x) iog(^) 



r<°)= 7 + ^(-l/2)+ ^ V2 \ 
/j 1} -1/2) 



a 



J 

! " -( i ^.9^2 J/2 r(-l/2)/(- 1) (-l/2)), (25) 



Here is a constant of the dimension of mass that is not fixed by the bare theory. It is worth mentioning 
that can be absorbed by the dimensional constant \i via its rescaling: ft = /ie r<,) / 2 . Physically, the cutoff 
at fi originates from inapplicability of our model at low energies due to trigonal warping and/or many-body 
effects as was discussed above. 

2. 5. Evaluation of the ultraviolet divergent terms 

Below we describe how the zeta - regularized effective action appears in a more conventional expression 
for the fermionic determinant. Namely, we consider 

- L eff = -^rTrlog' 



2T * V (id T ) 2 + Hi 
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^g s g v J™ dtt*- 1 (Tre-^ 2 +^ - Tre^f^ ^ 

r°° / o 7 r°° 



-v 2 r J t 



8vT0F 



L ' *=g,<ir I <nr —yiti^e ' •— • • 1 jjj- / <lr, • • • )|.._„ 

n>J 

UL jLg s9v v{2B) 1+J ' 2 T(s - 1/2) Fj(s - l/2)| s=0) 



Fj(s) = — / dtt s - X (^ e-^r* + J/2 - / dre- Jf ) (20) 
l s J Jo n>J Jo 



Here Hq is the Hamiltonian of the system in the absence of external magnetic field. It is worth mentioning 
that in this expression we do not omit zero modes of the Hamiltonians % and "Ho as these zero modes are 
to cancel each other. 

At J = 1 Eq. (l26t is convergent, and we find numerically using MAPLE package that F{— 1/2) = £(—1/2). 
This means that for the case J = 1 zeta - regularized effective action appears when the contribution of the 
free fermions is subtracted. It is worth mentioning that in the resulting zeta - regularized expression the zero 
modes of T~L are omitted while in Eq. (f2l)f are not. This means that the zero mode contribution is subtracted 
automatically when we subtract the contribution of free fermions. 

For the cases J > 2 Eq. (|26|) contains divergences: 

-L eff = -^—r=g a 9v v(2B) 1+J / 2 T(s — 1/2) Fj(s — 1/2, <5)| s= o, 

1 POO pOO 

Fj(s,S) = — / dtf- 1 ^ e-ra 4 + J/2- / dre~ T *) = fj(s) + F?*(s, 5) (27) 

Here 5 — - , where A is the ultraviolet cutoff. The divergent terms Fj lv (s, S) are calculated in Appendix 
D. It occurs that the zeta regularized effective action at J > 3 appears as a subdominant term. The values 
of the effective action for J = 1, 10 (together with the corrections due to the small electric field) are 
presented in Table [TJ For J = 2 the numerical factor in front of the term ~ vB 2 log ^ given in Eq. (6), Eq. 
(8) of 0] is reproduced with the magnetic scale A = V2vB, see also Ref. jj^]. For J > 3 we keep only the 
dominant ultraviolet divergent terms proportional to B 2 and the finite correction due to the electric field. 

The constants in front of the ultraviolet divergent terms at J > 3 do not have much sense because these 
terms appear to be of the same order as the terms that come from the large values of energy (> A), where 
the energy depends on the momentum (in the absence of external fields) as £ ~ \p\. We do not evaluate the 
latter terms here. According to P, [ljl[2(| we have v = A 1_J Vp. Therefore, the leading terms for J > 3 are 
~ v F B 2 A~ 1 . The complete effective action with the subdominant non - analytical terms and the divergent 
terms is given by Eq. (27) with the values of fj(s) and F^ m (s, 5) is presented in Appendix C and Appendix 
D. 



3. Effective action in the presence of magnetic field and small electric field 

3.1. Schrddinger equation 

Let us now consider the case when weak in-plane electric field (along the x-axis) is added. We consider 
the external vector potential in the form: A x = Et + By. The one-particle Hamiltonian in a subsequent 
parametrization has the form (cf. Eq.([2])) 

I ((p x + Et + By)-ip y 

H = v j 

y i^(p x + Et + By) + ipy) 

Again, we try the wave function as ^(t, x)) = ^ e %v * v i\) v {t). 
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3.2. Zero order approximation and Hall conductivity 
Similar to subsection 12.21 we change the variables: 



and denote 



A = F E - , u = VB(y+^±^) (29) 
a f = [~d u + u] , a = —= [d u + u] (30) 



Instead of Eq. ([6]) the Schrodinger equation reads 

i[d t + vB J / 2 \(a-a+)}ij 1 = v{2B) J ' 2 [a + ] J ^ 2 

i[d t + vB J ' 2 \{a- a + )]ip 2 = v(2B) J / 2 [a] J ^ 1 (31) 
When E is small, in the zeroth order approximation we have the analogue of Eq. ([7]): 

idt^x = v(2B) J ' 2 [a)] J ^ 2 

id t i> 2 = v(2B) J ' 2 [a] J ^ 1 (32) 

We have the same Landau levels as without electric field. However, the centers of orbits y c = — [p x + Et]/B 
move now slowly with the velocity v c = E/B along the y-axis. This means that the Hall effect takes place, 
i.e. when electric field along x axis is turned on, the current along y-axis appears if some of the energy states 
are occupied. The conductivity at zero temperature can be calculated as 

" B E nr.. c \ 9s9v (J 



9*9v l^^-B-56{n-£ n ) = —[-+I{}i) 



EL x L y ysyv ^ 2nE B ^ " ; 2tt 



Here fi > is the chemical potential. It has been taken into account that the zero energy levels are half-filled 
at (i — 0. The number of nonzero Landau levels with £ n < /i is denoted by I(^). We come to the well-known 
conclusion (see Ref.[l[) that at \l = +0 the Hall conductivity is equal to 



9 s 9v 



J_ 

4tt 



(33) 



3.3. The first order correction 

In the next approximation we take into account the second term in the l.h.s. of Eq. ([3T 

[e + i\(a-a))]^ 1 =2 J l 2 [a)] J ^ 2 
[e + i\(a- at)]V>2 =2 J / 2 [a] J ip 1 



(34) 



i\(a 



At - 



2 J ' 2 \a 



a) 2 J / 2 [a)] J 
J iA(at - a) 

Here £ = vB J / 2 e is the one - particle energy in the reference frame moving with the velocity E/B. In 
zero order approximation its value is given in subsection 12.21 

One can easily see that the first-order term in expansion of e in powers of A vanishes for n > J. For 
n < J — 1 there is the first order correction to the energy for J > 2. This correction is given by the eigenvalues 
xj^A of the J x J matrix i\Q with 



n 



( 


1 







-1 


V2 




... 






... 


-y/j-2 


V 


... 














(35) 
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In particular, we have: 

x (2) =±i; _ 

X (3) =0,±V3 

X (4) = ±2.334414217, ±0.7419637843 

X (5) = 0, ±2.856970013, ±1.355626179 

X (6) = ±3.324257431, ±0.6167065895, ±1.889175877 

X (7) = 0, ±3.750439718, ±1.154405395, ±2.366759409 

X (8) = ±4.144547191, ±0.5390798102, ±1.636519041, ±2.802485863 

X (9) = 0, ±4.512745858, ±1.023255666, ±2.076847980, ±3.205429001 

x (io) = ±4.859462833, ±3.581823478, ±0.4849357082, ±1.465989092, ±2.484325840 (36) 

3.4- The second order approximation 

One can easily see that the first-order term in expansion of e in powers of A vanishes for n > J and the 
second-order term 



e n = e (°) + A 2 £ (0) 1 , M @ - a)M 2 (37) 



should be considered. Here we use that at |n| > J 



/ B 1/4 „-u z /2g /, 

^ - V2 { S ign(n) ^g^ e-^/ 2ffw _ j( ,) ) ™ 



whereas at \n\ < J 



ip n =( 2l4^ e u2 ' 2R U u ) ) (39) 







For n> J Eq. (|3"T1) can be rewritten as 

A 2 / 1 



Ae„ = Cn -eW = — ( (Q) [Vn + 1 + - J + l] 2 + (Q) (Q) [v^+Vn-J] s 

e « — e n+l e ™ — £ n-l 

„„ ^ 1 (5T-[^ + T - + (Q \ 0) lVn--Vn—Jf) = A 2 ry r! 

£ ™ + £ (n+l) £ ™ +£ (n-l) 

(J-4)(2n-(J-l)) 



2J(n)i /2 2 J / 2 
In particular, in the case J = 1 we have 



(40) 



3A 2 



For n < J the second order correction vanishes. 

In order to make the consideration of the one - particle spectrum complete, we consider in Appendix B 
the problem in the gauge, where the electric field is introduced via the scalar potential Aq = Ex. In this 
gauge there are real energy levels £ n ,p y = jjPyi where £ n are calculated above while p y is the momentum 
along the y axis. The Galilean transformation of energy gives the values of the energy levels equal to £ n in 
the reference frame moving along the axis orthogonal to E with the velocity E/B. 
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3.5. Effective lagrangian 

In principle, in the presence of external field E the one - particle problem is not stationary, and we do 
not have usual energy levels in the original reference frame. However, we do have such levels in the frame 
moving with the velocity E/B along the axis orthogonal to the direction of electric field. In this case the 
expression for the effective lagrangian is derived in Appendix A. It occurs that this expression coincides with 
the usual one ^2 |£|/2, where the summation is over the energy levels of the system defined in the moving 
reference frame. This allows us to calculate the effective lagrangian in this case. 

For J = 1 the correction to the effective lagrangian reads: 

- A £ .„ . ± . \^Lf„ (41, 

where we consider Yl^Li as C( — 1/2) ■ 

For J > 2 the dominant contribution is given by the term linear in E due to the splitting of the lowest 
Landau Level: 

- - sgh^.B"' 't^m^i^ = -^.*^ 1/2 £ ix.i m 

7i—0 n—0 

For J > 4 the second order correction can be calculated without use of any regularization: 

Ar(2) BL x L y n j/ 2 ^ E 2 L x L y T E 2 ~ (J - 4)(2n - (J - 1)) 

_ eff = -^9.wB'i g^^pr = —^9**-zm g 2J(n) vv/ a 

L x L y E 2 (J-4)(2 gJ (-l/2)-(J-l)/j(l/2)) 

= ~^r^^m ^jj- 2 (43) 



Here in addition to the function fj given by Eq. (|22|) we use the function gj. These functions are the 
generalizations of the Riemann zeta-function and are given by: 

n>.J [U) -J n>J W-J 

(n)_j = n(n-l)...(n- J+l) (44) 

For J = 4 the correction in Eq. f4"3")l vanishes. For J > 4 expressions of Eq. f4"4"|) are convergent. 

In order to calculate the second order correction to the effective action at J — 1, 2, 3 we need to apply a 
certain regularization. First, let us consider the zeta - regularized expression. It can be calculated as follows: 



n 47r V^ vr ( s ) Jo 



\s=0 



v u n>J 

L , L 



y= gs g v Bd s (^-v^(2 B y^ 



2ra(J-4) + (5 J - J 2 — A) 



7 4i ^ 2 - BJ+1 



2 J +!J 

n>.J 



\s=0 



(45) 



This is the generalization of Eq. (|20p to the case of nonzero electric field. The sum over n is convergent due 
to the exponential factor. Coefficients r\ n arc defined in Eq. (|55"|) . We obtain {L e ff — L x L y l e ff); 

-1® = vB^ J / 2 +^ d s (^v- 2 %2By Js ^g s g v 2 J / 2 T(s- 1/2) f r ( S - 1/2)), (46) 



11 



-Al (2) - 1 E 



ds (fT7) v ' 2s{2Br Js ^ht 9s9v2J/2r{s + 1/2) 



eff ~ » WJ^ S \W) ' ' 

/ (J-1KJ-4) 2(J-4) xx 

I — /j( ' s + 1/2) ~ ^Trr- 9j(s ~ 1/2) J J 

Here the first line reproduces Eq. (12T1) as it should for the case of vanishing electric field. 

Similar to the case of pure magnetic field it is necessary first to calculate fj(s) and gj(s) for the values 
of s, where these series are convergent. Finally we must continue analytically the obtained functions of s 
to s = —1/2. In analogy to the theory of zeta-functions [25, .26], we shall find integral representations for 
these functions that are convergent at all values of s (see Appendix C). Then the resulting integrals can be 
evaluated numerically In addition to representation Eq. (|23[) we have the similar one for the function gj 
(s^O): 

fj (a ±1/2) a f -" ){±1/2) +ff\±l/2), 9J (s ± 1/2) « ^ (±1/2) + gf (±1/2) (47) 
s s 

The calculation of the values of fj 1: °\±l/2), g\ 1: °\±l/2) is described in Appendix C. It has been 

shown that fj ^(±1/2) may be nonzero for J = 2(21f + 1),K € Z while gj 1/2) may be nonzero for 

J = 2,4. We also have /j~ 1} (l/2) = for J > 3. Remarkably, for J = 2 the contributions of / 2 (_1) (l/2) and 

g[j 1 (— 1/2) to the logarithmic term (~ logB) cancel each other. Therefore, we have l e ff = l^)f + A/ e // 
with 

-Al ff - ^ t (0) 

7J - (^ 9 , 9 ,^/'r(-l/2)( ' J - ^ 4 » /',°»(l/2) - ji^i) b f(-l/2)l)), (48) 

Our results for the effective action at 1 < J < 10 are accumulated in Tabled] (These results are obtained 
using the data from Table [5]). From this table it follows that the ~ E 2 correction due to electric field is 
decreased fast with the increase of the number of layers. At the same time, the coefficient before the main 
term due to the magnetic field is increased. It is worth mentioning that at J = 4 the ~ E 2 correction to the 
effective action vanishes in this approximation. The linear term in electric field describes the linear Stark 
effect, which leads to the spontaneous (broken symmetry) electric polarization. This term can be presented 
as a scalar product E ■ P, where the vector P is directed along the spontaneous polarization, which in the 
presence of electric field is oriented along the field. 

3.6. Conventional regularization 

According to our experience due to the consideration of the case when there is the magnetic field only, the 
ultraviolet divergent terms may be present in the effective action even if its zeta - regularized version is finite. 
Therefore, let us consider the conventional regularization for the suspicious cases J = 1,2,3. We subtract 
the contribution of the fermions at B = 0, E = 0. However, approaching to this limit is performed along the 
line of constant A. This means that the boundary conditions for the fermionic fields in the functional integral 
for the free fermions are anti - periodic in the reference frame moving with the velocity V = Xv2 1 / 2 B^ J - 1 ^ 2 , 
(B — > 0) along the axis orthogonal to E. (Remind that for the system in the presence of external fields E 
and B the antiperiodic boundary conditions are adopted in the reference frame moving with the velocity 
E/B.) 



"eff 



9s9v 
' 2T 



L1 .'o 



Trlog ( 



(id r ) 



[id T ) 



H 



^±g s g vV (2B)^/ 2 



dtt s - 3 / 2 



n>J 



expl — (n)-jt — tX 



a (J -4)(2n -(■/-!)) - 
2 J J 
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n>J 



ex P (-(n)-jt)-J™dTex P (-T J t~t\ 2 {J 4) %/ J + f°° dr exp(-r J i)) | s=0 

(49) 



This is the generalization of Eq. (|26|) to the case of nonzero E. Here T-Lq is the Hamiltonian of the system 
at B, E — > but with the same A = -0172^7+1772 as in %. 

At J = 1 Eq. ((33) is convergent and is equal to the zeta regularized expression. This means that for 
the case J = 1 the zeta - regularized effective action appears when the contribution of the fermions with 
vanishing B and E is subtracted. The limit E, B — > is obtained along the line of constant A. 

At J = 2 there is the logarithmic ultraviolet divergence in Eq. (|4"9")l . At J = 3 Eq. (j4T)|) contains the 
divergency but only in L^^. 

3. 7. Magnetoelectic effect 

Let us remind (see Appendix A) that 

- L eff [E, B] = F[E, B] = -L x L y l ef f (50) 

where F[E, B] is the free energy of the system calculated in the reference frame, where electric field vanishes. 
It is equal to the effective lagrangian with the minus sign. The quantity 

F = F + L x L y h(B 2 — E 2 )/2 (51) 

(where h is the thickness of the graphene sheet) can be considered as the thermodynamical potential for 
fixed E and B. Actually, F is the effective lagrangian (with the minus sign) for the system of graphene and 
the constant electromagnetic field. As usual, we may introduce vectors D and H : 

D = P + E, H = B — M, (52) 

where P and M are electric and magnetic polarizations (of a unit volume) : 

dF 

-j—j—r = -PdE - MdB, (53) 
L x Lyh 

Then we have: 

— l —rdF = -DdE + HdB (54) 
L x Lyh 

Thermodynamical potential U with respect to variables B and D is related to F as follows: 

U = F + DE (55) 

Its differential is 

dU = EdD + HdB (56) 
L x Lyh 

One can see that this is related to the lagrangian L e ff = — F in the same way as the classical Hamiltonian 
(that is another definition of energy for the electromagnetic field): 

U = F ' E % = Ed ^E~~ leSh L eff = Leff- L x L y h(B 2 -E 2 )/2 (57) 

For the case of weak electric field the thermodynamical potential F per unit area depends on E and B 

as 2 

j^j-F[E, B] = j^-F[E = 0, B] - ujB^E - Ij-^jn (58) 



where the coefficients are given in Table [T] 
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The magnetization M = — j~L h fg anc ^ electric polarization P = — L X L h ^ of undoped multilayer 
graphene at zero temperature can be found from this expression. Due to the effect of electric field on 
Landau levels, the magnctoclcctric effect arises, that is, the dependence of magnetization on the electric field 
and electric polarization on the magnetic field. It is characterized by the quantity 

BM BP 1 d 2 F (1 „_ n/ , T E 

h—— = h 



BE BB L x L y BEBB 



It follows from Table [T] that the second term in the magnetoelectric effect becomes weaker when the 
number of layers is increased while the first one is increased with the increase of J. 

4. The system in the presence of electric field and small magnetic field 

4-1- Semiclassics in one-particle Schrodinger equation 

To establish the relation with the previous work [7| let us consider the gauge At = —Ex, A y = Bx. Then 
we have stationary Schrodinger equation = with 




(p x - i(p y + Bxfj 



(60) 



Ex 

We proceed with the rescaling z = I j x 7 u> — (7^7) e, 7 = B( J , and II = 

1 

("§) J+ Pv- Then instead of Eq. ([31]) we have: 

(z - u)i>i + {-iB z - i(n + iz)y ] i> 2 = 

{z-u)^ 2 + {-iB z +i{U + 1 z)) J i)) l =0 (61) 
The first-order semiclassical approximation for tpi^ gives 

$ = e ±l I (-(n+7^) 2 +(2-") 2/J ) 1/2 ^ ( 62 ) 

We denote here u = (z — uj) 2 / ,] /8 2 and 6 = II + juj. Then 

4 = e T&' +1 ij((i+ie J - 1 u'"->f-uf>-u->/^du (63) 

Integration over the classically forbidden region (l+70 J ~ 1 it J//2 ) 2 > u gives usthe pair production probability. 
In the limit 7 = (B = 0) we have 



hoi 2 



ao = 2JB(-,-) (64) 

The first-order perturbation in B results in [3] : 

■ 12 „- Qo e J+1 -7 2 Qie 3J - 1 



\V 

a 1 =3J 2 B(iH) (65 ) 
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4-2. Field-theoretical consideration 

The fact that in our approximation the particles do not interact with each other allows to reduce the 
field-theoretical problem to the quantum-mechanical one. Namely, we arrive at the following pattern. Modes 
for different values of momenta propagate independently. At t < to all states with negative values of energy 
are occupied while all states with positive values of energy are vacant. Their evolution in time is governed by 
the one-particle Schrodinger equation. At t = tg + T the wave function already has the nonzero component 
corresponding to positive energy. Its squared absolute value is the probability that the electron - hole pair 
is created. In this section we imply that the gauge is chosen such that A x — Et, A y = Bx. However, the 
probability that the electron-hole pair is created at the definite value of the momentum p y does not depend 
on the gauge chosen. Therefore we can use here the results of the previous subsection. 

Let us calculate the probability that the vacuum remains vacuum P v (vacuum persistence probability). 
According to the above presented calculation this probability is 



Pv = 

Px ,Py 



n(i-i ? , + i 2 ) 9s9 w 2ims (66) 



Here S is the effective action, the factors g s — 2 and g v — 2 are spin and valley degeneracies. The product 
is over the momenta that satisfy 

ET/2 >p x > -ET/2 (67) 
The total probability of the pair creation per unit area per unit time is 

2ImS" E 



= -9s9v- 



TL T L V 2ttL 

PV— Ly 



E lo g(i-b/+l 2 ) 



-9*£- I ^log(l-h + | 2 ) 
2n \v/ ^-^ n I 2tt 



The final result reads 



/1\7TT E^ s( J + 2 \r( 1 ^ 

9s9v \v) 2(J+l)7r 2 (a ) 1/(J+1) VJ+l/ li 

1 "" ° c(7$f) r (jTT) 



a = 2Js(|,^), ai =3J 2 i?(i,3J/2) (69) 

According to Ref. [31( a different quantity is considered as the pair production rate: 
T = (\v+\ 2 )/(L x L y T) 



J + 2 



1\jtt E^ / 1 w d2 ^A 4/(j+1) "i r (^n" 



2(J+l)7r 2 (a ) 1/( ' /+1) VJ+lVV \E) a { a J ~ 1)l[J+1) 



«r- 1)l{J+1) r( 7 i I ) 



5. Conclusions 



In this paper we calculated for the first time the effective Euler-Heisenberg action for the multilayer 
graphene at ABC stacking in the presence of external electric and magnetic fields at E « B. In the 
opposite limit B « E we calculated only the imaginary part of the effective action. The considered 
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effective field model is a kind of the quantum field theory with the anisotropic scaling r — > br, t — > b J t 
that is now becoming relevant. The particular anisotropic scaling with J = 3 has been applied by Hofava 
for construction of the quantum theory of gravity, which does not suffer from the ultraviolet divergences. 
Graphene and graphene like materials may serve as the condensed matter realization of the anisotropic 
scaling with arbitrary J. These materials have nodes in the fermionic spectrum, which are characterized by 
the integer momentum-space topological invariant N. Close to such a node fermions behave as 2+1 massless 
Dirac particles with energy spectrum £(p) = ±vp N . These fermions induce the terms in the action for 
electrodynamic fields, which obey the anisotropic scaling with J = N . 

Here we considered the fixed space dimension D — 2, but with point node of arbitrary topological charge 
N. This is somehow orthogonal to the relativistic systems studied in literature, which in case of massless 
fermions corresponds to the fixed J = 1, but with arbitrary space dimension D. In our case some features 
(in the presence of the external fields E « B) look similar, but some are new: 

1. For D = 2 and J = 2 the logarithmic term appears that is naturally expected from the one- loop consid- 
eration. In relativistic theories the same logarithmic term naturally appears in one-loop consideration 
in conventional electrodynamics with massless fermions, which corresponds to J = 1 and D = 3. 

2. For J > 3 there are terms that are divergent stronger than logarithmically. These terms depend on 
magnetic field but do not depend on electric field. 

3. At J > 2 the magnetoelectric effect is dominated by the lowest Landau level. Its degeneracy in 
the absence of electric field is J. In the presence of electric field the degeneracy is eliminated. The 
corresponding term in the effective action is proportional to E, which produces the analog of the Stark 
effect. The lowest subdominant terms are proportional to E 2 and are decreased fast with the increase 
of J. There is a specific value J = 4, at which the quadratic E 2 term is absent. 

4. The term which contains the logarithm ~ 8t InB appears only for special values of J, such as 
J = 2, 6, 10. This situation is similar to what occurs in relativistic theories, where the existence of 
the term ~ B~^~ In B depends on the space dimension D 11]. It would be interesting to consider the 
general case of arbitrary D and J. 

It is worth mentioning that at J < 3 the zeta regularization of the effective action at E « B gives the 
correct result. At the same time, for J > 3 the zeta regularization gives only the subdominant terms. The 
dominant ultraviolet divergent terms are calculated as well. This case is similar to the case of Dirac fermions 
in the presence of external fields E « B and with the dimension of space D > 3. 

In the case E >> B we calculate the imaginary part of the effective action with the small correction 
~B 2 . 

In future it will be instructive to consider electrodynamics arising in general case of arbitrary space 
dimension D in the vicinity of the manifold of zeroes in the fermionic energy spectrum of different dimensions. 
It appears that the point nodes in D = 3 described by topological charge N > 1 gives rise to the more 
complicated structure of the induced electromagnetic action: the QED has different scaling laws for different 
directions in space. For example, fermions emerging near the Weyl point with topological charge N = 2 
have the J = 1 scaling law for s pec trum along an anisotropy axis and the J = 2 scaling for the transverse 
directions (see Sec. 12.4 in Ref. |16j). 
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Appendix A: Effective action and one-particle spectrum 

In expression Eq. ([1]) for the effective action H is the one - particle hamiltonian in the presence of external 
electric and magnetic fields E and B. This hamiltonian may depend on time explicitly (as, for example, in 
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Eq. (|28p). T —> oo is time. It is implied that E is small enough, so the electron - hole pairs are not created. 
Our main supposition here is that there exists the transformation y — > y = y + Vt, with some velocity V 
such that in new variables (x,y,t) we have [idt — H[E, B,t]]ip(x,y,t) = [idt — H[E, B]]ip(x,y,i), where H 
does not depend on time. For Eq. (|28|) this is achieved for V = E/B. 

The system is considered with anti - periodic in time boundary conditions (in this new coordinates): 
ip(t + T,x,y) — —ip{t,x,y). Suppose we find the solution Q of the equation (id t — H)C = such that 
(s(t + T) = e~ l£T C,g- In this case £ is the eigenvalue of H and the analogue of the energy level. Actually, it 
is the energy level in the case, when T-L does not depend on time from the very beginning, and V = 0. Then 
^k,£ = e l T( 2fc+1 )*+ l£ *( £ is the eigenfunction of the operator (id t — H): 

{id t - n)^k, £ = -( J(2* + 1) + £)^ k ,s (7i) 



The product over k can be calculated as follows 



r7T 



iifc=_tf,...,tf[-(2fc + i) + £] = nfc=_ w ,„., w -(2* + i)) n fc [i 



ST 



Irp \— , ~, • J "T V ' J y ' ^(2fc+l) J 

« (^) 2Ar (nf = _ A , +1 (2/c + l))cos(£T/2), (72) 

where T = 2iVa, and a is the lattice spacing. This results in 

Det (iSt —H) = e- n °n n cos(£„T/2), (73) 

where f2o depends on the details of the regularization but does not depend neither on T nor on the spectrum 
in continuum limit. The values £ n depend on the parameters of the hamiltonian, index n enumerates these 
values. 

Partition function receives the form (see also [l3, 37|): 

Z(T) = e- n ° J2 exp(i^J2 £ ^- iT J2 Kn£n ) = e ~"° J2 exp(-iT^ K n £ n ) (74) 

{K„}=0,1 n n {Jf„}=0,l n 



Following [13( we interpret Eq. (|74l) as follows. represents the number of occupied states with the 
values of "energy" £ n , K n may be 0, 1. The term J2 n £ n vanishes because the values £ n come in pairs with 
opposite signs. 

After the Wick rotation we arrive at 

Z(-it) = e~ n ° J2 exp(-r^X„£„) (75) 

{if„}=0,l n 

Here r = where T is temperature. In the formal limit T — > only X n = [1 — sign£ n ]/2 survives. 
Thus we get 



r 

n,£ n <0 



Z(-i/T) = e-° exp(-i £ n[E,B]), T^O (76) 

Comparing Eq d76t and Eq. (p} we obtain 



Z(T) = e- no exp(-zT £ B]) (77) 

n,£„<0 

and 

L e// [£,B] = - ^ £ n [-E, £] = +(1/2) ^ \£n[E, B]\ = —F[E, B] (78) 

n,£ n <0 n, 

Here F[E, B] is the free energy of the system in the presence of constant external fields E and B in the 
reference frame moving with the velocity E/B in the direction orthogonal to the direction of E. (Actually, 
in this reference frame the electric field is absent.) S — TL e ff[E, B] is the effective action. 
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Appendix B: Energy levels in the presence of magnetic field and small electric field 



5.1. Schrodinger equation 

Here we consider the case when weak in-plane electric field (along the x-axis) is added. We use here the 
same gauge as in Sect. 14.11 Namely, the external vector potential has the form: A y — Bx while the scalar 
potential is A t = Ex. The one-particle Hamiltonian has the form of Eq. (1ST))) 



H 



Ex 



v(p x + i{p v + Bx)j 
We try the wave function as ty(t, x)) = J2 Py e lPvV ^ p (x,t). 



>(px - i{Py + Bx ) 
Ex 



(79) 



5.2. Zero order approximation 
Let us change the variables: 



£ + §Py 

£ — T In 1 ^ 



£ 



vB J/2 ' 



vB J / 2 



A = 



E 



V2vB( J + 1 )/ 2 ' 



= VB(x + ^), u c = 



B 



Py 



and denote 



at 



-j=[-d u +u], a=-^=[d u + u] 



(80) 
(81) 



Instead of Eq. (f3"Tj) the Schrodinger equation reads 
eip = (e — \V2u c )ip = 



A(at+fl) (-i) J 2 J / 2 [a] J 



When E is small, in the zeroth order approximation over A we have: 

(£ + §P„)Vi - v{2B) J '*[tf] J ^ 

(E + ^ Pv )tP 2 = v{2B) J / 2 [a] J ^ 1 
is 



(82) 



(We cannot neglect here ^p v because p y may be large.) We have the same Landau levels as without electric 
field shifted by -§p y . 

5.3. The first order correction 

In the next approximation we have 



A(at + a) (-») J 2 J / 2 [at] J 
2 ^ J A(at + a) 



i 3 2 J l-\a 



4> 



One can easily see that the first-order term in expansion of e in powers of A vanishes for n > J. For 
n < J — 1 there is the first order correction to the energy for J > 2. This correction is given by the eigenvalues 



(J) 



A of the J x J matrix \tt with 



fi = 



/ 1 

1 V2 

... 

o ... V7^2 o 

V o ... o VJ~\ 



o 



(83) 



We have the same eigenvalues as in Eq. (|36l) . 
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5.4- The second order approximation 

One can easily see that the first-order term in expansion of e in powers of A vanishes for n > J and the 
second-order term 

e n = ?£» + A 2 £ _ (Q) 1 _ (0) \^ + a)M 2 (84) 

should be considered. Here we use that the functions ip n are given by Eq. , Eq. (fTUj) . 
For n> J Eq. (fMf can be rewritten as 

Ae„ = A 2 7?„, 

(J-4)(2n-(J-l)) 
2J(n)i / J V/ 2 

For n < J the second order correction vanishes. 

Appendix C: Calculation of the functions fj(s) and gj(s). 

5.5. Integral representations 

In this section we describe regular procedure for the calculation of /j(s) and gj(s). 
For s > 1/J (s > 2/ J) we have 

M*) = £[(^3p "to = E [(n) _" ]a+ i (86) 



We apply to these sums the Plana summation formula (2||, Vol. 1, 1.9(11): 

°° i />oo />oo 

E/W = o/(°)+/ f(r)dr + i (f(it)-f(-it) 



n=0 

This formula works if: 



1. /(a;) is regular for Rex > 0, 

2. e" 27r|t| /(r + it) -> at t -> oo, < t < oo, 
3- / e _27r l t l|/(r + it)\dt — > at r — > oo 



Therefore, we obtain: 

1 /"" /' v / \ <// 



(87) 



M») = ^7j]7+/ o (r + J)Z S jdT + ij ((it + J) = --(-it + J) = ^- 5 - F _ I) S >l/J 

OO 



sjOO = 2[7^ + y o (T + J):r 1 (^ + ^)^ 



Here 



+*/ o (^ + J):fr 1 (^ + i)-(-^ + ^:r 1 Hi + i))^r- T; *>2/J (88) 



(r + J)_, 7 = E^+Z+l) =(r + J )( r + j_i)... (r + i) (89 ) 
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Now the divergences of the sums over n at s < 1/J and s < 2/ J correspondingly are concentrated within 
the integrals over r. At s >1/J(s>2/J) we may represent these integrals as follows: 

/>oo />oo 

Jj(s) = / (r + J)!} dr, Ij(s) = (r + J)Z s ,-\t + J) dr 
Jo Jo 

= j> zZjdz = 2>(«)(1 + e- 17rJs ) + y zZjdz 

0= <j> zZj^zdz = Zj(s)(l - e -^J(«+i)) + (J- l) e - i7rJ ( s+1 )/ J (s + 1) + J zZj^zdz (90) 

Here contour C consists of the integral from — oo + iO to — 1 + iO, the part of the circle that belongs to 
the upper half - plane, starts at — 1 + iO, and ends at J + iO, then the part of the real axis from J + iO to 
oo + iO. The contour is closed via the half - circle at infinity in the upper half - plane. We obtain: 

i r J 

7 ->( S ) = - 1 + e -i*Js f Z -J dZ 



-1 



Ij(s) 



1 



1 — e _i7rJ ( s + 1 ) / i 



j 



--J 



( J _ l ) p -iirJ(s+l) 



Here integrals § are taken along the half - circle z = + ^-^-e luJ , ui € [0, n]. The given integrals are 
convergent for all values of s, probably, except for s — (2K + 1)/J and s = 2Q/J — 1, K, Q G Z. 

For the ordinary Riemann zeta function there is the Hermit representation ( 26], Vol.1, 1.10, (7)). Using 
the expressions listed above we derive the analogue of this representation: 

M _ _L 1 f J f°° sin s [arctg -| + ... + arctg t] dt 

Ij[S) ~ 2[J!] S l + e~ i « J 'f_ 1 Z - J J (t 2 + J 2 )*/ 2 ...(t 2 + l)*/ 2 e 2 -*-l' 

J 1 f J , J-l e -*rJ(«+i) r- 7 

ffj(s) = 2[J!F~ l-e-^ +1 ) l/- 7 + 1 = e-^W) i + f_ x Z -J dz 

f°° sin( (g + l)[arctg j, + ... + arctg t] - arctg fl) (f 2 + J 2 ) 1 / 2 g 

7o (* 2 + J 2 )( s+1 )/ 2 -.(i 2 + l)( s+1 )/ 2 e 27rt -l 1 j 

This expression gives the analytical continuation of the series Eq. (f86|) to all values of s. From this repre- 
sentation we conclude that fj(s) may have a simple pole at s = (2K + 1)/J, K € Z while gj(s) may have 
poles at s = (2iT + 1)/J or s = (2Q + 1)/J - 1. 

5.#. Particular cases 

At J = 1 expression (ji?2"j) is the Hermit integral representation for the £ - function (see (2||, Vol.1, 1.10, 
(7))- 

We have (s 0): 

/j(s ±l/2) « /r ) (±l/2) +/ ( 0)(±1/2); gj(a±1/2) ^ ^" 1) (±V2) +g (0) (±1/2) (93) 

We have verified integral representations Eq. (|92p using MAPLE package for several values of s and J such 
that the series Eq. (|86[) are convergent. Next, we calculated numerically the values of ff' ^(±1/2) and 
9j {— 1/2) using these integral representations at J = 2, 10. 



At J = 2 we come to: 



1 , T r~\ , n f°° sin s [ arc *g f + arctgi] d< 



/2(S) ~ 2 s +i +/2(s) + 2 io (t 2 + 2 2 ) s / 2 (i 2 + l)-/ 2 e 2 - t -l' 



1 - f s f°° sin((s + 1) [arctg | + arctgi] -arctg |]) (t 2 + 2 2 ) 1 l 2 dt 



*(*) - V +HS) + 2 J ( t2 + 22 )( s+1 1/ 2 (, 2 + 1)(s+1 )/ 2 e ^_! (9 4 ) 



20 



For J = 2 we may check expression Eq. (|91|) using integral representation for the hypergeometric function 
([il, (2.12), Eq. (5)): 

/2(S) = 2^( g -l/2) F(S ' X ' ^ V2) = 2^(.s-l/2) + Ts\ + ^TTT) i + 
Therefore, we calculate: 



sin(- (1/2) (arctg (t/2) + arctgt)) dt 



/•OO 

-2 



(4 + *a)-i/4(l + t 2 )-V4 e 27rt 
= -0.3321609172 

^<-i» - -^(i + fi + fi + ..., = -o.oe 25 , 

/f^) = ^+d s [^F(.s, 1,2s, 1/2)] | s=1/2 

sin((l/2)(arctg (t/2) + arctgt)) dt 



OO 



(4 + t 2 ) 1 /4(l +i 2)l/4 g 27 r 

= 0.01902854064, 

= 2^^(1/2,1,1, 1/2) = 1/2, 

= ^(1/2,1, 1,1/2) = 1/4, 
^(-5) = V2 + a s [ £ 2 |^F( S ,l,2 S + l,l/2)]| s= _ 1/2 



-2 



sin((l/2)(-arctg(t/2) + arctgt)) dt 



(4 + t3)-V*(l+t3)V4 e 2jrt -l 
= -0.8676318824, (96) 

For J > 2 we proceed in a different way. For example, for the case J — 3 we have 

11 f 00 sins[arctg 4 + arctg | + arctgt] 



o (t 2 + 3 2 ) s / 2 (t 2 + 2 2 ) s / 2 (t 2 + l) s / 2 e 27rt - 1 ' 
13 - r sin( (a + 1) [arctg | + arctg § + arctgt] - arctg |]) ft 2 + 3 2 ) 1 / 2 

53l.s; O + 3 W + /. , 2 o 2 \( s +l)/2("*2 , Q2\( s +l)/2^2 , -{\<s+l)/2 e 2nt _ 1 ^ ^ 



2 6;5 +i J w y Q (t2 + 32)(«+i)/a( t a +2 2)(»+i)/2( t a + 1 )(«+i)/2 
Here ^(s) and ^(s) are to be calculated using Eq. (19"T1) . We get: 

/J 0) (-~) = 0.1580218749, / 3 (0) (^) = 1.687512381, gi°\-^) = -0.576826000, (98) 
At J = 6, 10 the residues f(~V{-l/2) are calculated as follows: 

f^i-m = -±-j> J x V*dz (99) 

f (o) 



The values /} u; (-l/2) are calculated as 



sin (-1/2) [arctg j + ... + arctgt] dt 



^ ' (t 2 + J 2 )- 1 /4...( i 2 + 1) -i/4 e^rt.j ( 10 °) 
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J 








/r j (D 


^(-§) 




1 


-0.2078862251 





-1.460354509 





-0.2078862251 





2 


-0.3321609172 


-0.0625 


0.01902854064 


0.5 


-0.8676318824 


0.25 


3 


0.1580218749 





1.687512381 





-0.576826000 





4 


-0.2742668544 





0.5508597242 





0.3077772742 


0.25 


5 


-1.650991958 





0.2018069599 





1.792805053 





6 


-1.299466810 


-0.2460937500 


0.07309406118 





0.5888090809 





7 


1.789507452 





0.02547240488 





0.2149754417 





8 


-10.74171499 





0.008486905628 





0.07746594906 





9 


-73.39264066 





0.002702793021 





0.02686060136 





10 


-28.2691703 


-7.922363278 


0.0008241640781 





0.008909168673 






Table 2: The values of /j(±l/2) and gj(-l/2). 



In a similar way we have calculated the other values collected in Table [5J For J > 4 the original 
expression for gj(— 1/2) Eq. (jgrJl is convergent, and we use it to obtain the corresponding values. Therefore, 
5^ 1} (-l/2) = 0, J > 4. In a similar way /j(l/2) is convergent for J > 2 and /j~ 1} (l/2) = 0, J > 2. At the 

same time it is not excluded that fj 1 ' (—1/2) ^ at J = 6, 10, 14, .... Indeed, we calculate /j" 1 ^— 1/2) for 
J = 6, 10 and find that these values do not vanish. 



Appendix D: Calculation of the ultraviolet divergent terms in F( — 1/2). 

5. 7. The case J = 2 
For J = 2 we have: 



1 A 00 / ' Z" 00 \ 

F2M) = f^)/- d " e"^ T * + 2/2-y o dre^ 2 ') (101) 



with <5 = -^nr-, where A is the ultraviolet cutoff. We rewrite this expression using Plana's summation 



formula as follows 



1 Z" 00 /I Z" 00 Z" 00 

_J_ / dtt--ifi[ e -a*-l]+ / d T e-( T+1 )( r+2 )* + 3/2- / dre- 



+2 



sin s [arctg | + arctg t] dt 



r( 



o (t 2 + 2 2 ) s / 2 {t 2 + l) s / 2 e 2wt -1 



1 , „ f°° sin s [arctg | + arctg t] 



2(2!) s J (t 2 + 2 2 y/ 2 (t 2 + l) s / 2 e 2 ^ -1 : 

At s = —1/2 the second integral over t is convergent while in the first one there is the divergent term 
that corresponds to the integration over the region of small t: 

Ft{-l/%6) « i ^^^'''dtt- 1 / 2 - 1 (it 1 /2r(l/ 2 )) Rj llog«5 (103) 

This expression gives the logarithmic term in the effective action after the renormalization of effective charge 
(that results in the change A — > fi). At the same time at — 1/2 < s < expression Eq. (|102[) is convergent, 
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J 


Ff v (-1/2,5) 


1 


— 


2 




3 


~ TTTi 

-0.7639792361 <T 1/6 


4 


-0.8642091734 <T 1/4 


5 


-1.094743796 <r 3/1 ° 


6 


-1.393109122 6~ 1/a + 0.2460937500 log 8 


7 


-1.746814339 <T 5/12 - 6.155346286 <T 1/12 


8 


-2.151696783 cT 3/8 - 5.842427394 8~ 1/8 


9 


-2.605972745 x <r 7/18 - 6.875813125 x <T 1/6 


10 


-3.108750471 x <T 2/5 - 8.590800550 x cT 1/5 + 7.922363278 log 5 



Table 3: Divergent terms in Fj(— 1/2,5). Here <5 = - ^ 2 f^ , where A is the ultraviolet cutoff. 



but has the pole at s = — 1/2 with the residue —1/16. In this region we have: 

r0+ 



1 



T(s)[l - , 
1 



2ttis] 



dtt s 



-i (3/2 + dr [e -^)(r + W) _ e -^) 



2(2!) s 



sin s [arctg | + arctg t) dt 
o (^ + 2 2 ) s / 2 (i 2 + l) s /2 e 2 "*-!' 



(104) 



Here integral f is over the closed path that starts at infinity, goes to zero along the real axis, then encloses 
t = clockwise and comes back to infinity along the real axis. This contour can be deformed in such a 
way that it is placed at infinity and is wrapped around t — clockwise. Written in this form Eq. (|104l) is 
convergent for all values of s except for s = 0, ±1, ±2, .... For s > 1/2 we come to expression of Eq. (1M1) and 
thus derive F 2 (s,0) = f 2 (s),s > -1/2 . 



5.8. The case J = 3 
For J = 3 we have: 



F 3 (s,6) 



T(s) 



dtt s 



e T^aTT* + 3/2 - 



n>2 



(105) 



with 5 = — A a '- , where A is the ultraviolet cutoff. We rewrite this expression using Plana's summation 
formula as follows: 



F 3 (s,6) 



i 1 poo 

dT e -(r +t V3 )(T+2t V3 )(r+ 3 tl / 3) + ^ _ 1 / ^ ^ 

t ' JO 



1 



sins [arctg 4 + arctg | + arctg t] 



2(3!) s J Q (i 2 + 3 2 ) s / 2 (i 2 + 2 2 ) s / 2 (t 2 + l) s / 2 e 27rt -l' 
The second integral over t is convergent while in the first one there is the divergent term 

'v 2 (2B) 3 \-i/6 



Ft" (-1/2,6) = 



r(- 



/3 



0.7639792361 

1/6 



(106) 



(107) 



This gives the divergent term in the effective action ~ v (2_B) 1+3 / 2 ( - ^ ~ _B 2 A 1 



/3^2/3^ 



5.9. Arbitrary J 

For arbitrary J we have: 



r( s ) 



■(E 



e + J/2 - I dre~ 

o 



(108) 
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v (2B) 

with 5 — — A2 ; , where A is the ultraviolet cutoff. We rewrite this expression using Plana's summation 
formula as follows: 

1 r sins[arctg j + ... + arctg£] dt 

2(J!) S i (*2 + J 2 ) s / 2 ...(t 2 + l) s / 2 e 27r *-l' 1 J 

The second integral over t is convergent while in the hrst one there may appear the divergent terms 

Fp{-i/2,s) = —^ r d tt-v^( Yi A '< tK/j ) 

^ ' ' JS l<K<J/2 

« J2 A K S~^ K ^+A J/2 \og5 (110) 

\<K<J/2 

Here the log term may exist only for even J. The divergent terms for 1 < J < 10 are represented in Table [3] 
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